Stroboscopic theory of atomic statistics in the micromaser 
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We study the statistics of the atoms emerging from the cavity of a micromaser in a dynamical, 
discrete-time 'stroboscopic' description which takes into account the measurements made, in general, 
with imperfect efficiencies t; < 1, on the states of the outcoming atoms. Inverted atoms enter 
stochastically, in general, with a binomial distribution in discrete time; but we also consider the 
continuous-time limit of this input statistics which is Poissonian. We envisage two alternative 
experimental procedures: one of these is to consider a fixed number A'^ of atoms pumped into the 
cavity and subsequently leaving it to undergo state detection; the other is to consider input of the 
excited atoms and their subsequent detection and collection in a fixed time t. We consider, in 
particular, the steady state behaviors achieved in the two limits, N ^ oo and t —> oo, as well as 
the approaches to these two limits. Although these limits are the same for the state of the cavity 
field, they are not the same, in general, for the observable outcoming atom statistics. We evaluate, 
in particular, Mandel's Q-parameters Qe (Qg) for outcoming atoms detected in their excited states 
(ground states), for both N -^ oo and i ^ oo, as functions of N^x = RTc'. R is the mean rate 
of entry for the incoming atoms and Tc is the cavity damping time. The behavior of these atomic 
Q-parameters is compared with that parameter for the cavity field. 
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I. INTRODUCTION 

The micromaser, in which two-level atoms enter a high- 
er cavity one at a time and interact with a single mode 
of the cavity field before subsequently leaving the cav- 
ity, provides a valuable means of testing various aspects 
of the quantized-field-quantized-atom interaction [Ij-H ; 
while the one-atom-one-mode model of Jaynes and Cum- 
mings [|l5| used to describe the strictly one-atom micro- 
maser, as well as the A^-atom-one-mode model |1^ for 
those cases where more than one atom is in the cavity at 
the same time, are of fundamental interest in the theory 
of exactly solvable quantum models |ll| , |l7| -[l9| . Despite, 
for example, a recent operator equations of motions study 
of the dynamics of these micromasers [|9|, the body of 
theoretical investigations has been largely based on the 
characteristic properties of the cavity field averaged over 
the states of the atoms leaving the cavity. However, be- 
cause of the practical difficulties associated with measur- 
ing the states of this cavity field, these theoretical predic- 
tions can not be checked experimentally directly. In the 
experiments, the measurements are actually made on the 
states of the atoms leaving the cavity [HI 14 20-02]. The 
properties of the cavity field must then be inferred from 
the statistics of the states of the outcoming atoms. It is, 
therefore, essential to know the relationship between the 
field statistics and this atomic statistics. This relation- 
ship, evidently, requires the knowledge of not only the 
field statistics but also that of the atomic statistics. 

The atomic statistics in the micromaser, namely the 
outcoming atom statistics, refers to the probability dis- 
tribution functions for observing the atoms in one state 
or the other after they have left the cavity. These func- 
tions can obviously be found by repeatedly monitoring 



these states, either for a fixed number A^ of the atoms 
pumped into the cavity and subsequently leaving it or 
for the atoms pumped for a fixed duration t of time and 
subsequently collected in that time. The two descriptions 
are evidently equivalent if the time interval between the 
successive atoms entering the cavity is a fixed number but 
not if that time interval is a random variable. The dy- 
namics of the cavity field can similarly be described either 
as a function of the number A^ of atoms traversing the 
cavity or as a function of total passage time t. Of partic- 
ular interest is the asymptotic, i.e., the steady state, be- 
haviors of the system. The meaning of these asymptotic 
limits must, of course, depend upon the chosen descrip- 
tion: the limit corresponds to the limit A^ — )■ oo if the 
dynamics is described in terms of the number of atoms 
and it corresponds to the limit i — )■ oo if the dynamics 
is described in terms of the collection time. However, it 
has already been shown |H| that the density matrix of the 
cavity field in both of these limits is the same. 

The atomic statistics of the micromaser, on the other 
hand, has only been studied in terms of the atoms col- 
lected in a fixed interval of time and that too in the 
'coarse-grained' description of the equations of motion 
for the cavity field density operator pj [ni 23 - 29] . How- 
ever, and despite the fact that coarse-graining appears to 
involve a coarse averaging over the successive entries of 
the atoms [^ 23 , it has been known ]Q that when the 
probability distribution of the intervals between succes- 
sive atoms entering the cavity is exponential, i.e., when 
the pumping mechanism is Poissonian, coarse-graining 
which assumes a steady rate R (say) of entry for the 
atoms, is precisely equivalent to assuming this strictly 
Poissonian pumping of the atoms with the same mean 
rate R. The refs. |29-E9] for example adopt this coarse- 



grained description from the outset and then collect the 
outcoming atoms in a fixed interval of time. It is there- 
fore of interest to see whether or not these two alterna- 
tive descriptions of a fixed number N of atoms collected 
or of a fixed time t for that collection are asymptoti- 
cally equivalent in terms of the outcoming atom statis- 
tics rather than the field statistics for a general pumping 
statistics which includes Poissonian pumping only as a 
special case. Moreover, and despite the fundamental in- 
terest of this problem to ergodic theory, insofar as these 
two limits prove to be different (and they do) there is the 
possibility at the level of actual experiment of gaining ad- 
ditional information by collecting at fixed N on the one 
hand and at fixed t on the other. We shall see that this 
additional information can be helpful in determining the 
actual state of the cavity field in the micromaser. Thus 
it is this general question, of fixed N or of fixed t, which 
is the question addressed in this paper. 

In this paper we find the probability of the occupation 
of the two atomic states (to be called |e), excited and 
\g), ground) for a fixed number N of atoms leaving the 
cavity and we find the same probability for the group 
of atoms collected in a fixed time t. We shall assume 
the pumping mechanism is binomial. Both regular and 
Poissonian pumping are then special cases of this bino- 
mial pumping. We shall evaluate the parameters Q for 
these two different cases: Q takes the form of Mandel's 
'Q-parameter' [Q and is a measure of the deviation of 
the variance for the population of a given state from the 
value it would have if the distribution of that population 
were Poissonian. We shall find that these Q-parameters 
for each of the two states are generally different in the 
two asymptotic limits, the limits N —^ oo and t —^ oo. 
We investigate the dependence of these differences on the 
parameter N^x where Ng^ is the average number of atoms 
entering the cavity in one decay time of that cavity. The 
Q-parameters Q^ for the excited state |e) are found to 
be different for all Ng^ in these two asymptotic limits. 
However, the Q-parameters Qg for the lower state \g) are 
different in the two asymptotic limits for small N^x but 
those differences gradually disappear with an increase in 
Nex- The value of N^x after which the differences be- 
tween the two asymptotic results become insignificant 
are found to depend upon the interaction time. At zero 
temperature it is Qe for a fixed number N of atoms as 
N ^ oo which closely follows Qf the Q-parameter for the 
cavity field: but at finite temperatures this close relation- 
ship changes significantly. These observations concern 
only the low lying trapping states: for other equilibrium 
states of the cavity field Qe joins substantially in its be- 
havior with the behavior of the other Q's even at zero 
temperature. 

The paper is organized as follows:- By following the 
stroboscopic description, we derive in section || expres- 
sions for the joint probability of finding sets of succes- 
sively entering atoms in a particular sequence of detected 
states. We also derive a number of alternative strobo- 
scopic maps connecting successive states of the cavity 



field, and which have taken into account as an ensem- 
ble average the different sequences of atomic state de- 
tections. From these in the continuous time i.e., Poisson 
limit of the input statistics we derive master equations for 
the field density operator for fixed N. From this result 
we can regain the usual coarse-grained master equation 
for the field density operator. In the section HI we use 
the joint probability found in section || to find the joint 
probability of detecting given numbers of atoms in each 
of the two states on the passage of N active atoms and 
from it the variances in the number of atoms detected 
in a given state. We also find the joint probability of 
detecting given numbers of atoms in the two states for 
a fixed collection time t of the atoms and from this de- 
termine the variances for the atoms collected in the fixed 
time t. Numerical results comparing the two descriptions 
are then presented in the section IV . In an Appendix we 
outline a method for calculating Q^ and Qg for fixed N 
and the corresponding Qe and Qg for fixed t. 



II. JOINT PROBABILITY FOR STATE 
DETECTION OF SUCCESSIVE ATOMS 

In this section we derive an expression for the joint 
probability for detecting each of a set of successively en- 
tering atoms in either one of their two states after they 
have left the micromaser cavity. 

We consider the usual micromaser system in which a 
high-Q cavity is pumped by Rydberg atoms at a rate 
R so low that almost always there is at most one atom 
at a time in the cavity. We assume that the atoms enter 
the cavity prepared in the state |e) which is coupled reso- 
nantly by a cavity mode of frequency wq to a lower energy 
state \g) which, for the sake of convenience, we refer to 
as the ground state. The interaction between such an 
atom and the e.m. field as the atom traverses the cavity 
is governed by the Jaynes-Cummings Hamiltonian Pq| 



H 



iLJoa*a + LUaSz + g[a S- + S+aj , (2.1) 



where {a, a'') are the field annihilation and creation op- 
erators which obey the commutation relation [a, a^] — 
1 for bosons (they satisfy the Heisenberg-Weyl alge- 
bra |l|,0); operators S+ = \e){g\,S- = \g){e\,S, = 
(l/2)[|e)(e| — |(7)((7|] are the atomic operators (and sat- 
isfy the su(2) Lie algebra for total spin S ^ ^ [ |ll| , |l9| ), 
and g is the atom-field coupling constant. The cavity 
field and the atoms interact also with a heat-bath held 
at a constant temperature. This coupling induces the 
excited atomic state to decay spontaneously and to the 
decay of the cavity field. However, the time for the spon- 
taneous decay between two Rydberg levels, which are 
the levels of interest here, is long enough compared with 
the convenient choices of the times of fiight tint of the 
atoms through the cavity to enable one to ignore, to a 
very good approximation, the effects of the atomic cou- 
pling to the thermal reservoir {tint is the interaction time 



between the atoms and the cavity field and is ~ 35/x s 
IPOJ). Moreover the couphng of the cavity field to the 
heat-bath may also be ignored during this atom-transit 
time tint because this time is several orders of magni- 
tude smaller than the decay time Tc of the cavity field 
(Tc ^ 0.2 s [Q). Together these two acceptable approx- 
imations mean that the dynamical evolutio n is simply 
unitary under the Hamiltonian H given by (^j) during 
the short time tint- Note that the micromaser system 
has been studied numerically and in depth during this 
short atomic transit time tint with the coupling of the 
field to the heat-bath included |lC-|lJ] and it has also 
been studied analytically |Tlp^ - 34 1 . A naly tical results 
for the Jaynes-Cummings Hamiltonian (l2.1 | ) c oupled to 
the heat-bath have also been given in ref. |35[| in a dif- 
ferent context. All of these studies together confirm the 
fact that the coupling to the heat-bath does not play 
any significant role if tint is short enough compared with 
the atomic and field damping times. As noted this is 
the case, for example, in the experiments of ref. pfl] — 
although ||ll[ reports some small differences (of the order 
of percent) from Meystre's formula pi for the probabil- 
ity of finding n photons in the cavity mode even when 
N = 6000 atoms have passed through the cavity for pa- 
rameters within this experimental range. The numerical 
work ||ll| is essentially exact and the differences of the 
order of percent are compatible with the errors involved 
in the second of the two approximations itemized above 
together with the approximation ti^i — ti — tint ~ ti+i — ti 
made below; however N = 6000 atoms may still not be 
close enough to equilibrium. 

Still, whatever the exact situation for the short time 
tint , the coupling of the field with the heat-bath must cer- 
tainly never be ignored during the time intervals between 
the arrivals of two consecutive atoms into the cavity — for 
these time intervals are usually significant as compared 
with the damping time of the field. If the heat-bath has 
an average of n thermal photons at the cavity frequency 
luq and if 2k = T^^ is the rate of loss of photons then 
the evolution of the density matrix pf of the field due 
to its interaction with the thermal reservoir is governed 
by the usual master equation taken in the rotating frame 
(cf. e.g., |lj] and references) which is 



dPf 
dt 



K[{n + l){2apfa' — a'apj — pfa'a) 

+n{2a' pfG — aa^ Pf — pfaa')] 
Lpf. 



(2.2) 



Note that L is norm preserving i.e., -^Trpf ~ TrLpf = 
0, for the right side is expressible as a sum of commuta- 
tors the characteristic form for Markovian dynamics |36| ] 
(actually the sum in (2.2) is a sum of double commuta- 
tors). Now, let Pf{ti) be the density matrix of the cavity 
field at the time ti of the entry of the i*'* atom into the 
cavity. Since, in this paper, the atoms are assumed to en- 
ter the cavity in their excited state |e), it follows that the 
state of the combined system of the atom and the field at 



the time ti is the outer product of pf{ti) and |e)(e|. The 
system then evolves, as discussed above, under the action 
of the Hamiltonian (2.1) alone for the time tint — the time 
of fiight of the atom through the cavity. It can then be 
shown that under this unitary evolution during t^t the 
density matrix p(ti + tint) of the combined system as the 



■th 



atom leaves the cavity is given by 



Piti+tmt) = Fe{tint)\e){e\ + Fg{tint)\g){g\ 

+F,git,nt)\e){g\+Fl{t„,t)\9){e\\pf{t,) 

(2.3) 

where 



Fe{tint)pfiti) = cos{gtintVa^a + l)pf{ti) 
X cos{gtintVa^a + 1), 

Fg{tint)Pf{ti) = a^ 



Vfltfl+l 



sm 



{9tintVa^a + l)pf(ti) 



X sva.{gtint\/ a^ a + 1) 



1 



Vatfl -t- 1 



a. 



Feg{tint)Pf{ti) = icos{gt,ntVa^a + l)pf{ti) 

X sin(ff^;;„t\/atQ + l) a. (2.4) 

VaTa + 1 



(2.5) 



Significant quantities are, e.g. |l|, |ll| , |23| , |24p2| , p5[ 
(3m = sin^ (gtintVm) = 1 - a™. 



The quantity f3m+i is the a priori probability that an ex- 
cited atom after traversing the cavity containing m pho- 
tons would exit it in its ground state and am+i is the 
probability of exiting in its excited state (cf. e.g., |p7|). 
The expression (2.3) shows that the actual probability 
P{e;ti) {P{g;ti)) that the i*^ atom exits the cavity in 
the excited (ground) state is 



P(e;t,)=Tr/[FeP/(iO] 
(p(g;tO=Tr/[Fgp;(t,)] 



(2.6) 



where Trf denotes the operation of trace over the field. 
After leaving the cavity, the atom passes through a de- 
tector which determines its state. We assume that the 
time that the atom takes to arrive at the detector after 
leaving the cavity is small enough for any losses that the 
atom may suffer due to spontaneous emission can be ig- 
nored. However, the probability of detection of the atom 
in either state may still not be the same as the probabil- 
ity with which it leaves the cavity in that state. This is 
because the detector may not be of unit efficiency. Fol- 
lowing ref. ||21|], for example, we assume that the detector 
has an efficiency rje {rjg) for detecting an atom in its ex- 
cited (ground) state. Clearly, the probability Pd{e;ti) 
{Pd{g;ti)) that the i*'' atom is detected in the state |e) 
{\g)) after leaving the cavity is then given by 



= Tr[F,^,dPf{U)], v^=e,g (2.7) 



where P(e; ti) {P{g; ti)) is the probabihty, given by (p.6[), 
that the atom leaves the cavity in the state |e) {\g)) and 



^ u.d ^U^ U 7 



e,5- 



(2.8) 



Here, we have introduced the subscript d on those prob- 
abihties and operators which depend on the d etec tor effi- 
ciency and F^^d, like F^, depends on tint as in (2.4). Since 
the detector efficiency is not necessarily unity, there is a 
finite probability that the atom passes through the detec- 
tor without its state being detected. It is straightforward 
to see that the probability Pd{n;ti) that the atom goes 
through without any state detection is 



where 



Pd{n;U) = Tr[Fn^dPf{U)], 



Fn,d = il-Ve)Fe + il-Va)Fg, 



(2.9) 



(2.10) 



and depends on tint- We have thus evaluated the proba- 
bility of detec tion in a given state by (2.7), or of no state 
detection by (2^), for an atom as it comes out of the 
cavity. These expressions (2/7) and (^^) show that the 
probability of detecting an atom in one of its states or 
that of no detection is obtained by operating on the den- 
sity matrix (density operator) of the field, taken at the 
time of the entry of that atom, by the operator F^^d (now 
with 1/ = e,g,n) followed by the operation of trace over 
the field. Next we use this procedure to find the prob- 
ability of a particular outcome of the process of state 
detection for the next, i.e. the {i + 1)*'*, atom. 

Now, in order to determine the probability of detection 
or of no detection of a state of the (z -I- 1)*'' atom on its 
exit from the cavity, we need to know the state of field 
at the time i^+i of its entry into the cavity. Since the i*'' 
atom is subjected to detection, the state of the field in 
the cavity at the time of entry of the {i + 1)*'' atom is 
conveniently described by the conditional density opera- 
tor pf{ti+i\i'i; ti) |p8[ which characterizes the field at the 
time ti+i of entry of the {i + 1)*'' atom into the cavity 
under the condition that the i*'' atom, which entered at 
ti, is detected at ti + tint in the excited state Ui — e, 
is detected in the ground state Vi — g, or goes through 
without any state detection Vi = n. By following the pro- 
cedure outhned following equation (2.10) the probability 



of a particular outcome for the process of detection on 
the {i -\- 1)*'' atom is then found to be 

Pd{vi+i\ti+i\vi;ti) ^Trj F^^^^^dPfiU+iWhU) ■ (2.11) 

Next we must find the expression for the conditional den- 
sity operator p/(ii_|_i|t'i;ii) in terms of the state of the 
field at the time ti of the entry of the i*^ atom. This 
expression is found by first determining the conditional 



density operator pf{ti-\-tint\vi] ti) of the field at the time 
ti + tint of the exit of the z"* atom. 



By comparison with (2.3) it is straightforward to see 
that the state of the cavity field at the time of exit of the 
i*'' atom corresponding to each outcome of its detection 
is (with dependence on tint included) 



Pf{tt +t^nt\vi,ti) 



Fvi^d{tint)Pf{ti) 
Pd{Vi\ti) 



Vi = e,g,n. 



(2.12) 



After the exit of the i atom, the density operator 



( p.l2[ ) evolves according to the master equation (2.2) 
for the time interval i^+i — ti — tint until the time t^+i 
when the (i -I- l)*'* atom enters. So far in all quan- 
titative micromaser experiments one chooses as far as 
is possible |3^ that i^+i — i^ » tint- Hence we let 
ti^i—ti — tint ~ ti^i—ti = Ti. In this way the state of the 
cavity field at the time i^+i of the entry of the (i -I- 1)*'* 
atom is given in terms of the state at the earlier time 
ti + tint and in terms of the outcome of the measurement 
on the atom i by 

Pf{ti+iWi;ti) ^ D{n)pf{ti + tintWi;ti), (2.13) 

where 



Dixi) ^ exp(LTi), 



(2.14) 



is the formal solution of ( pT^). On substituting pf{ti 
tintWi',ti) from (2.12) into (p. 13) we thus get 



, I . DiT,)F^^^d{t,nt)pfit,) 

Pfit.^MA) = j,^^^;^^ ,-. 



e,g,n; 



(2.15) 



for the density operator of the field at the time of the 
entry of the {i + 1)*'* atom under the condition that the 
i*'* atom exits the cavity in the state ji^i), Vi = e,g or goes 
undetected (vi = n). 

The Eq. ( 2.15| ) shows that the state of the field at the 
time of entry of the {i + 1)*'' atom corresponding to a 
particular outcome of the process of measurement on the 
i*'' atom is given by operating on the density matrix at 
the time of entry of the i*^ atom by the operator F^-^d 
[vi = e,g,?i) followed by operation with the operator 
D and then normalizing the resulting expression to unit 
trace. 

The desired expression for the conditional probability 
Pd{i'i+i] ti+i\i>i; ti) of atomic det ection may now be ob- 
tained by substituting ( 2.15| ) in (2.11). This expression 
therefore reads 



Pd{vi+l'-,ti+i\vi]ti 



P'^f Fu,+x,dP>{n)Fu,4Pf{ti) 



Pd{i^i]ti) 



(2.16) 



In what follows we require the joint probability 
P(t'i+i; ii+i, Ui^ ti) that the outcome of the process of de- 
tection on the «*'' atom is Vi together with the outcome 



for the (i + 1)*'^ atom being fi+i. This probabiHty can 
be obtained by using the relationship 

Pd{.Vi+i\k+i,v.,\ti) = P{vi+i;ti+i\i'i]ti)Pd{i'i]ti), (2.17) 

between th e join t and the c onditional probabilities. On 
combining (2.16) and (2.17) we find the desired expres- 
sion 

Pd{vi+i]ti+i,Vi;ti) = Tr f[Fi,^^^^dD{Ti)Fu^^dPf{ti)]-, 

(2.18) 

for the joint probability of a particular pair of outcomes 
as a result of measurement on two successive atoms en- 
tering the cavity at times ti and i^+i. We can now repeat 
the preceding arguments and show that the joint prob- 
ability Pd{vN;tN, . . . ,iy2;t2,iyi;ti) = Pd({i^fc;tfe}Ar) that 
the outcome of measurements at their exits on the atoms 
entering the cavity at times ti,t2, ■ ■ ■ ,tj^ is vi,V2, ■ ■ ■ ,vn 
is given by 

Pd{{vk;tk}N) = Tr f[Fy^ ^dD{TN-i)Fut, -i,d ■ ■ ■ 

. . . D{T2)F,,^dD{Ti)F,,^dPfiti)l (2.19) 

where p/(ii) is the density operator of the field at the 
time of entry of first atom into the cavity. 

Sim ilarly, by following the procedure outlined after 
( 2.15 ), the state of the cavity field at the time t^ of the 
entry of the N*'^ atom is described by the conditional 
density operator 



Pf{tN\{l^k]tk}N-l) = 

[D(rAr_i)F^„_j^d. 



.D{T^)F,,^dPfiti)\ 



Pdi{i^k;tk}N- 



(2.20) 



Note that the density matrix Pf{tN) characterizing the 
state of the cavity field without any condition on the state 
of the exiting atoms must be related to the conditional 
state by the relation 

Pfi^N) = ^ Pf{tN\{l^k\tk}N~l)Pd{{vk;tk}N-l), 
{uk=e,g,n} 

(2.21) 

where the summation denotes the sum over all the pos- 
sible outcomes. On substituting (2.20|) in (^.21) we get 



{i/fc=e,g,n} 

.■■D{T2)F,,,dD{Ti)F,,^dPf{ti) 



(2.22) 



The summation on each of the Vk for one label k can now 
be carried out by noting that for each k 



E 



F, 



i^k,d 



F,_ 



F 



g.d 



F„ 



Fe + F„ 



(2.23) 



where the equality follows from the use of (2.7) and 
( p.lO|) . Note that this step at ( p. 23 ) neatly eliminates 
the detector efRciencies T^e, % from (2.22). For a first in- 
vestigation let us now assume that the atoms arrive with 



a regular spacing Tp, i.e., Ti is Tp for all i. The Eq. (2.22) 
can then be seen to be exactly 



Pf{tN) = [D{Tp)F^f-^pf{t, 



where 



Fn = Fe + F, 



(2.24) 



(2.25) 



Note from (2_^) that FQp{ti) is the density operator of 



the field at the time of exit of the i atom if the state 
of this outcoming atom is left undetermined. Moreover 
the Eq. (2.24) shows that, as is required by (2.22), the 



ensemble averaged field does not depend in any way upon 
the efficienci es of the detectors. 

The Eq. (2.24) is evidently equivalent to the strobo- 
scopic equation or 'map' 



PfiU 



D{Tp)FoPfit,), 



(2.26) 



for the regular input as it was used in e.g. ref. p3| . The 
route by which (2.26) is arrived at here is, of course, not 
at all the one that is usually followed to reach this ex- 
pression. The usual derivation (e.g. psl) docs not make 
any reference, from the outset, to the processes of atomic 
state measurement. The derivation given above specifi- 
cally demonstrates the consistency of that approach with 
that developed here now taking proper account of the 
process of atomic measurement. 

The derivation of (2.26) given above shows actually 
that, for regular inputs at least, all effects of the atomic 
state detection processes described by the efficiencies rji 
vanish under the ensemble average. Indeed because of 
(^.23|) these efficiencies vanish through the ensemble av- 
erage (2.22) for p/(ijv) for a^j/ statistics of the sequence r^ 
of successive atomic inputs. For the cavity field therefore, 
effects of state detection have vanished naturally from the 
ensemble averaged theory. Actual effects of the detec- 
tion process via the atoms on the state of the cavity field 
were envisaged already in [y and were investigated via 
numerical methods in |21 1; one of the present authors also 
investigated such effects numerically 11 40] by however 
studying only single realizations of the atomic inputs — 
where substantial effects can be discerned ^^l- Within 
coarse-grained theory (that is, equivalently, for continu- 
ous time Poisson inputs — see below) such effects of the 
atomic detection processes on the state of the cavity field 
were also investigated in [g5[ . Also within coarse-grained 
theory, by introducing a remarkable non-linear master 
equation taking account of the outcomes of the atomic 
measurements the ref. ^6| shows how these measurement 
process effects are ultimately eliminated via th e ens em- 
ble averaging within this theory. Our step at (2.23) on 
the ensemble averaged, and usual linear, density opera- 
tor for the field, shows that quite generally for discrete 



time entries of the atoms it is not possible to detect via 
the atoms any effect of the measurements on the atoms 
on the state of the cavity field. This last statement fol- 



lows from the fact that we show in section III how the 



variances for the observed atom statistics depend only 
on the ensemble averaged field density oper ator pf{tN)- 
Of course this expression Pf{tN), Eq. ( ^.22| ), must itself 
be ensemble averaged over the individual realizations of 
the sequ ence r^, of the atomic inputs. In the derivation 



of (2.24) we took the Ti to be constants, r^ = Tp for all 



I. We also assumed even in (2.22) that the interaction 
time tint was a constant. However, in actual experimen- 
tal situations, the time difference ti between the arrival 
of successive atoms will not be the same for all the atoms. 
And even the time of interaction tint of the atoms may 
vary because different atoms may travel with different 
speeds. In fact, in practice, the times Ti and tint rnust 
both be taken to be random variables. 

In this paper we shall consider only the variations in 
the pumping time by following the widely used model 
[0,|| of pumping time fluctuations. According to this 
model, atoms leave the source oven at regular intervals 
T (say) and are then prepared in a given excited state 
|e) with a certain probability p prior to their entering 
the cavity in that excited state. Thus atoms entering 
with regular spacings T have only the probability p of 
being in |e) and only atoms arriving in |e) are effective 
in determining the dynamics of the micromaser. Even 
so one may still describe this dynamics in terms of the 
number K of all the atoms or equivalently in terms of 
the discrete time tx = KT that all these K atoms take 
to arrive in the cavity. Alternatively the dynamics may 
be described in terms of the number N of atoms entering 
the cavity specifically in their active excited states |e). 

First we find the expression for the joint probability 
Pd{{^k] tk}N,tK) that TV atoms enter the cavity in their 
excited states |e) in time tx and that the outcome of the 
measurement on those atoms after they leave the cavity 
is i'i,V2, ■ . . jVn- To find this probability let us assume 
that the state of the cavity field at the time i = is Pf{0). 
Now, assume that the atoms numbered ki, i — 1,2, ..., N 
are excited to the respective state |e). The probabihty 
for this event is obtained by using the fact that the prob- 
ability of m failures followed by a success in a binomial 
process in which the probability of a success is p is given 
by p(l — p)™. The time interval between the arrival of 
active atoms numbered i — 1 and i (i = 1,2,..., N) is 
therefore given by {ki — ki-i)T with fco = 1. T he p rob- 
ability Pd{.{vk',tk\ n ,tK) can be obtained from ( 2.1£ ) (i) 
by replacing the damping operator D{Ti) between the 
atoms numbered i and i — 1 by D{[ki — ki-i)T), (ii) by 
noting that the field density operator Pf{ti) at the time 
of the entry of the first active atom is D{(ki — l)T)pf{0) 
and (iii) by performing the summation over all ki to ac- 
count for all possible realizations oi {ki}. It then follows 
that 



Pdi{i^k;tk}N,tK)^p^il-p) 



K fcjv-1 

kN=N kN-i=N-l 



ko-l 



..-Y^Trf \d{{K -kN + l)T)F,^,dD{{kN - kN-i)T) 
fci=i 

y-F,^_,^d . ■■D{{k2 - ki)T)F,,^dD{{ki - l)T)p}[Q) 

(2.27) 

where, for later manipulative convenience, we have multi- 
plied on the left by D{{K—kN+l)T) the factor appearing 
under the trace. This does not change the results because 
of the fact that L is norm preserving and Tr(Lp) = so 
that 



Tr{Dp) = Tr{p), 



(2.28) 



for any p. 

The evaluation of (2.27) is facilitated by defining the 
(probability) generating function 



f{N,x)^ Y, x''Pdi{iyk;tk}N,tK), 



K=N 



so that 

Pdi{vk;tk}N,tK) 



M;si^("->|- 



(2.29) 



(2.30) 



On u sing the expression (2.27) for Pd{{i^k',tk}N,tK) in 
( ^.29 ) and on carrying out the summations over the {ki} 
and K we get 

f{N,x) = {px)^Trf[Dpix)eMLT)F,^^dDpix)eMLT) 
xF,^_,,d . . . Dp{x)eMLT)F,,^dDpix)pf{0)], 

(2.31) 



where 



Dp{x) 



1- x(l-p)exp(LT)' 



(2.32) 



Note that the joint probability Pd{{i^k'-, ifclw) for the out- 
comes {vk} on the passage of N active atoms irrespective 
of the time tx, obtained from Pdi{i^k',tk}N,tK) by per- 
forming the summation over K, is given by 

oo 

Pd{{>yk;tk}N)^ Y, Pd{{'^k;tk}N,tK) = f{N,l), (2.33) 



where the last equality is by virtue of (2.29). The joint 
probability Pd{{i^k', tk}, tx) for the outcomes {I'k} in time 
tx irrespective of the number of active atoms is obtained, 
on the other hand, by performing the summation over N 
in Pd{{i^k;tk}N,tK), i.e., 

K 

Pd{{vk-M,tK) = Y Pd{{vk\tk}N,tK). (2.34) 

JV=0 



We will use ( 2.33| ) and (2.34) in the next section to find 
the probability that a given number of atoms are detected 
in a given state after the passage of a fixed number TV 
of active atoms irrespective of time or that for a fixed 
collection time t^- Here we demonstrate the consistency 
of the approach developed so far by deriving the equation 
of evolution for the cavity field averaged over all possible 
outcomes of detection. 

The cavity field averaged over all possible outcomes 
of detection at the time of entry of the if*'' atom when 
N atoms have entered the cavity in their active excited 
states |e) follows from (2.22) and the considerations out- 
lined above and can be seen to be given by 



{i^i— e,(/,n} kn — 1 kn-i — 1 

fe2-l 

. . . E [oiTiK + 1 - kM))F,^^dD{T{kM - kN-i)) 



fci=i 



xF,„_,,d . . . D{T{k2 - ki))F,,,dD{T[ki - l))p/(0) . 

(2.35) 



On ca rrying out the summation over Vi in (2.35) by using 
(^.23|) we get 



K *;„-! 



pHa^,^k)=p^(i-p)^'-^E E ••■ 

fe2-l 

E [d{T(K + 1 - kN))FoD{T{kN - kN-ij) 



ki=i 



xFo . . . DiTik2 - ki))FoD{T{ki - l))p/(0) 



(2.36) 



where Fo is again as defined in ( 2.25 ) and the averaging 
over the Ui has again elimina ted t he effects of the atomic 
detection processes. From ( 2.36 ) it is now straightfor- 
ward to derive the recurrence relation 



Pf{N,tK+i) = exp(Lr) il-p)pfiN,tK) 



+pFoPf{N~l,tK) 



(2.37) 



which in principle determines the ensemble averaged cav- 
ity field pf(N,tK) at the time tx under the condition 
that in that time N atoms have entered the cavity in 
their active excited states le). Note once more that this 



ensemble averaged cavity field, relation (2.37) is indepen- 
dent of the detector efficiencies. It can be seen to reduce 



to ( 2.26 ) for reg ular inputs by setting p — I and T — Tp. 
The relation ( 2.37 ) can now be used to determine the 
field density matrix at the time t^ or at the time of entry 
of the A^*'' excited atom. The field density matrix Pf{tK) 
at the time tx is given by Pf{N, tx), summed over all N , 
whereas the field density matrix Pf{N) after the passage 



of a fixed number N of active atoms is given by pf(N,tK), 
summed o ver al l K. Now, on performing the summation 
over N in (2.37) we arrive exactly at the familiar map 



Pf{tK+i) = exp(LT) (1 - p) + pFo p}{tK), (2.38) 



already derived in ref . ||g] . The derivation of ( 2.38| ) in |^] 
was however quite different from that just given now; in 
g no account whatsoever was taken of the sets of mea- 
surement on the atoms; here these sets of measurement 
are fundamental to the whole analysis eventhough they 
play no role in the final result. In the same way the sum- 
mation over K on ( ^.37 ) together with a readjustment of 
the terms leads to the map 



Pf{N+l) 



peMLT)Fo 
l-(l-p)exp(LT) 



PfiN). 



(2.39) 



The relation ( ^.39 ) can also be derived from the map 
( p. 26 ) by assuming that the pumping time Tp in it is a 
random variable corresponding to the binomial process 
and on performing the average over all of its possible 
realizations. It is this procedure which was followed in 
ref. PI in order to derive the equation for the case of 
strictly Poisson pumping. 

Of particular inter est i s the steady state solution of 
t he ra ap s ([2.3^ ) and ( 2.39] ). The steady state of the map 
( ^.38|) (( p.39| )) is the state reached in the asymptotic limit 
Ik ^ oo (N -^ oo) and is such that pf(tK+i) = Pfi^K) 
{ pfjN + 1) = pf{N)). Hence the steady state p^s of 
(^.38|) is the solution of the operator equation 



1 - exp(LT) ({l-p)+ pF„ pss = 0. (2.40) 



Howe ver it is easy to see that the steady state of the map 
( p. 39 ) also obeys ( ^.40 ). Hence the steady state density 
matrix of the cavity field is the same whether the dynam- 
ics is described in terms of time or in terms of the number 
of active atoms traversing the cavity. The exact solution 
of ( 2.40| ) is known only for the case of Poisson pumping 
discussed next. 

A widely studied and experimentally important case 
of random pumping is Poisson pumping. This corre- 
sponds to the limit p — > 0, T — > of binomial pump- 
ing such that p/T -^ R where the constant R is the 
average rate of pumping of the active atoms. In this 
limit, the atoms enter the cavity continuously. In other 
words, in the Poisson limit, the discrete time tx be- 
comes a continuous variable t. Hence one can write 
Pf{N,tK+i) -pf{N,tK) = T{d/dt)pf{N,t). The recur- 
rence relation (2.37) then reduces to the equation 



j,PfiN,t) 



LpfiN,t) + R{Fo 



l)pfiN~l,t) , 

(2.41) 



coupling Pf{N, t) and pf{N — 1, t) and which on iteration 
determines the evolution of the field at time t under the 



condition that N active atoms enter the cavity in that 
time. The equation for the field density matrix Pf(t) at 
time t irrespec tive of the number of active atoms can be 
derived from ( 2.41 ) by performing the summation over 
N. This leads to the exact master equation 






L + RiFo-1) pf{t). 



(2.42) 



This master equation is the same as the one that is ob- 
tained in the coarse-grained approximation without mak- 
ing any reference to the pumping statistics (see for exam- 
ple |2^ or jn], [^ — wher e coa rse-grained averages are 
actually taken). The Eq. ( 2.42|) can also be derived by 
taking the Poisson limit of ( 2.3j ). Finally, the density 
matrix of the field after the passage of a fixed number 
N of active atoms, taken irrespective of the time for the 
case of Poissonian pumping, is found from ( 2.39| ) to be 
given by 



PfiN + l) = 



1-L/R 



FoPfiN). 



(2.43) 



This equation is the same as the one derived in ref. [|l| for 
t he P oisson pumping process achieved by starting from 
(2.26) and assuming the pumping times Tp to be ran- 
dom variables corresponding to a Poisson process and 
then performing the average over all realizations of this 

pumping time. 

Now cf. Eq. (2.27) and observe that we have there the 



exact stroboscopic expression for the probability that N 
atoms in their active excited state are pumped into the 
cavity in a given time tn and that the process of state 
detection after these atoms exit the cavity gives a partic- 
ular sequence of outcomes. We can use this expression 
to study various aspects of the atomic statistics includ- 
ing the correlations between the atoms exiting the cavity 
at different times. We shall study these correlations in 
a later paper. Here we restrict attention to the study of 
the probability of detecting a given number of atoms in 
a given state. 



III. DISTRIBUTION FUNCTION FOR ATOMIC 
STATE POPULATIONS 



In the last section we determined at Eq. (2.27) 
the joint probability Pd{{i^k', tk, }n, ^k) for the outcomes 
{vk = e,5,n} as a result of measurement on N excited 
atoms passing through the cavity in time Ik- Here we 
use that joint probability to determine the number distri- 
bution 'w{Ni.,Ng,N) which is the probability of detect- 
ing Nf, atoms in the excited state and Ng atoms in the 
ground state on the passage of A'' active atoms and also 
the number distribution wi^Ne, Ng,tK) for a fixed time 
interval tjf of observation. We use these expressions to 
find the variances in the number of atoms exiting the 
cavity in one state or the other. 



Recall that the probability Pd{{i^k',tk}N) is the prob- 
ability that the outcome of state measurement on N 
atoms passing through the cavity in an arbitrary time 
is t he set {vk} (yk = e.^g^n). This probability is given 
by ( ^.33D . Now, the probability w{Ne,Ng,N) that N^ 



atoms are detected in the state |e) and Ng in the state 
|g) on the passage of TV excited atoms is clearly the sum 
of all those joint probabilities Pd{{i^k', tk}N) {vk = e, g, n) 
in which N^ of the v^s correspond to the excited state 
and Ng of the v^s correspond to the ground state. We 



know from ( 2.33| ) th at -Pd di^fc; ^fcjjv) = f{N,l) where 
f{N,x) is given by ( |2.31 ). It is then straightforward 
to see that 'w{Ne,Ng,N) is the coefficient of y^'z^^i in 

the expansion of Trl LDp(l)pexp(LT)f F„_d -|- yFe^d + 
Dp{l)pf{0)> where Dp(x) is given by ( 2.32 ). 



zF, 



g,d 



In other words 



w{N,,Ng,N) 



^N.+N, 



-Tr 



N^\Ng\dy^-dz^o^'\ 
X \Fr,^d + yFe4 + zFg^d 



A„ 



N 



Pf 



(il)} 



z=0 

y=o' 
(3.1) 



where Pf{ti) = Dp{l)pf{0) is the density operator of the 
field at the time of entry of the first excited atom, 



Ap = Dp{l)pcxp{Lp/ R), 



(3.2) 



and we have substituted T = p/R where R is the average 
rate at whi ch th e atoms are excited. On substituting for 
Fn,d from ( 2.10 ) and on changing the var iabl es y and z 
to 1 — y and 1 — z respectively, one finds (3.1) becomes 



wiN,,Ng,N) 



(_)(We+Wg) J_N,+ 



Na 



NelNgl 



{Fo 



dy^'dz^'^ 
yr/eFe - zr]gFg 



Tri 



K. 



N 



Pf 



{h)} 



z=l 



(3.3) 



The distribution for the number of atoms detected in the 
ground state alone, obtained by summing ( |3.3| ) over N^, 
thus reads 



Wg{Ng,N) 



Ng\ dZ^f' 



Tr{[Ap(^Fo~zT]gFg 



N 



'Pf 



ih)} 



(3.4) 



together with a similar expression for We {N^ , N) in which 
the quantities with suffix g are replaced by the cor- 
responding ones with suffix e. Note that N^, and Ng 
are not independent variables, so their distributions are 
related. For the case of unit detection efficiency for 
each of the two states, N — N^ + Ng which implies 
We{Ne,N) = Wg{N — Ne,N). On using the easily verifi- 
able relation 



1 d'' 

1 



(a + xb] 



N 



x=0 



{N^ky.dx^-'' 



(xA + b\ 



N 



(3.5) 



it can then be verified that Wg{Ng,N) given by (3.4) 
and the corresponding expression for We{Ne,N) indeed 
satisfy the relation We{Ne,N) =Wg{N - Ne,N). 

Next we evaluate the atomic population distributions 
for the case when the observation is made for a fixed 
time Ik rather than for a fixed number of atoms. The 
joint probabihty Pd{{'^k',tk},tK) for the outcomes of de - 
tection in tha t cas e is determined by working from ( 2.34 ) 
rather than ( ^.33| ). The probability w{Ne,Ng,tK) that 
Ne atoms are detected in the state |e) and Ng in the state 
\g) in time Ik is clearly the sum over all those joint prob- 
abilities Pd{{vk]tk}, K) [vk = e,g,n) in which Ng of the 
ly's correspond to the excited state and Ng of the v's cor- 
respond to the ground state. Wit h Pd {ii^k-tk},tK) given 
by (p.34D taken with (|2.30|) and (|2.31|) it follows that 



w{Ne,Ng,tK) 



^N,+Ng+K 



NJNglKldy'^-dz'^odx^ 



E^P-)"" 



N 



Tr{ [Dp exp{Lp/R) (Fq - (1 - y)ri,Fe - (1 - z)TjgF, 
DpPfiO)} 



N 



x—y—z—0 



(3.6) 



Now, on performing the summation over N an d on rear- 
ranging the terms, one finds that the Eq. ( |3.6| ) yields 



1 



^N,+Ng + K 



w{N,,Ng,tK) - j^Jj^^lKldyN.dz^.dx^ 

xTr{ [l - xexp{Lp/R) (l - p + p{Fo - (1 - y)f]eFe 



(1 - ^)VgFg) 



Pf 



(0)} 



x—y—z—0 



(3.7) 



The differentiation with respect to x can now be per- 
formed trivially to get 



wiN,,Ng,t) = 



Tr 



NelNgl dyN-dz'^9 
[Fo - PivVeFe + zrjgFg 



A„ 



Rt/p 



Pf 



(0)} 



(3.8) 



where 



Ap = exp{Lp/R) , Fq = 1-p + pFo. 



(3.9) 



We have replaced ijc by t and substituted K — Ik /T — 
Rt/p. The probability of detecting Ng atoms in the 
ground state irrespective of the number of atoms in the 
upper state in time i, found by performing the summa- 
tion over Ne in (3.8), then reads 



Wg{Ng,t) 



Ngl dZ^9 



Tr 



a(Fo 



zpilgFg 



Rt/p 



'PfiO) 



(3.10) 



in which A = Ap and there is a similar expression for 
We ( -^e , t) obtained by replacing the quantities with the 
sufhx g by the corresponding ones with the suffix e. 

In the continuous time Poisson limit p — > we recover 
from (pT§) for fixed t the expression 



w{N„Ng,t) 



NelNgl dy^-dz 



rN, 



rr{exp[(L + i?(Fo-l) 



-yRrjeFe - zRrjgFg t 



Pf 



(0)} 
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(3.11) 



derived in ref. |26[| by their very different methods. The 
expressions for Wg{Ng,t) and We(iVe,t) given in |26l are 
similarly recovered in the Poisson limit of (|3.10|) and 
from the Poisson limit of the corresponding expression 

for We{Ne,t). 

Next we use the number distributions to find the vari- 
ances in the numbers of atoms detected in one state or 
the other. It is straightforward to see from the expres- 
sions for ■Wg{Ng, N) and We{Ne, t) that at fixed N for the 
ground states \g) 



{Ng)N = -j-Tr[[A(^Fo-zTjgFg 



N 



}pfitl) 



2=1 

(3.12) 



and that 

{Ng{Ng - \))^ 



d^ r r / M ^ 1 

—Tr[\t.[F^-ZT^gFg)\ p/(ti)| 



z=l 

(3.13) 



in which A = Ap while there are similar expressions for 
the excited states |e) obtained by replacing the quantities 
with the index g by the corresponding ones with the index 
e. The corresponding expressions for a fixed observation 
time t are given by 



Rtip 



{Ng)t = -^Tr[ [A(Fo - zprigFg)\ p/(0)} 



z=l 

(3.14) 



and 



{Ng{Ng - 1)), = ^Tr[ [a(Fo - zpr^gF, 



Rt/p 



P/(0)}| • 
(3.15) 



with si milar expressions for e rather than g. We use 
( |3.12 )-( |3.15|) to evaluate the parameters Qg defined in 
Mandel's ||31[] Q-parameter form by 



Qo 



{NgJNg - 1)) 

(Ng) 



(Ng). 



(3.16) 



Definition ( 3.1(j ) has a simil ar fo rm for Qe with e re- 
placing g. These quantities ( 3.16 ) are measures of the 
deviation of the observed distribution from a Poisson 
distribution. The values Qg — means that the cor- 
responding statistical distribution is Poissonian; Qg < 
and Qg > imply sub- and super-Poissonian distribu- 
tions respectively. An approximate relationship between 
Qg and Q f the corresponding Q-function for the cavity 
field both evaluated for fixed t is derived in ref . [Ell . How- 
ever this relation has since been found to be of limited 
validity 



lb 26 1 . The numerical calculations described in 
the next section for both fixed t and fixed N lead us to 
a similar conclusion. 



Now, the operations of differentiation in (3.12)-(3.1J;) 
can be performed by using the easily verifiable results 
that 



d_ 
dx 



(A + xB)" = ^(A + xBfB(A + a;5)"^'=-\ 



fc=0 



(3.17) 



and 



dx"^ 



[A + xBY = Y^[{A + xB)^B{A + xB] 



k-l-l 



xB{A + xB 

X B{A + xB 



k,l 



i'{A + xB)''B{A + xBy 
''-^-% (3.18) 



In this paper we shall evaluate (3.12)-( p.l5| ) by assuming, 
in accordance with ref. p6| , that the system is initially 
in the steady state , i.e. , p/(ii) = P/(0) = Pss where pss 
is the solution of (2.40). We make use of the fact that 
Dp and Fq are both trace conserving, i.e.. 



Trf[Dpp] = Trf[p] , Trf[Fop] - Trfip], 



(3.19) 



for an arbitrary p. On evaluating (3.12) and (3.13) and 
substituting the resulting expressions in (3.16) one finds 



{Ng)N^%NTrf[Fgp, 



(3.20) 



and 



2^a 



^' TrflFgp,,] 
1 



X 



1 - ApFo 



Trf[Fg[l~ 

^^p^ gPss 



l-(ApFo)^' 






N{l-K 
-VgNTrf[Fgp, 



(3.21) 



which is the desired expression for Qg with particular 
reference to the detection of atoms in the ground state 
\g) for a fixed number A'^ of excited atoms entering the 
cavity in an arbitrary time. Similarly, in the case when 
t he ob serva tion i s made for a fixed time t the expressions 
(|T4|) and ( ^l5| ) yield 



{N, 



git 



VgRiTrf[FgPs 



(3.22) 



and 



Qa^ 



^%P^ 



' Trf[Fgp. 
1 



,Trt 



/I l-(Apfo)^*/P \ 
'Vp i?i(l-A„Fo)^ 



X 



1 A fi j^pJ^gPss 
1-ApFo^ 



C(l iip-1 

TjgRtTrf[Fgp, 



(3.23) 



where the tilde on the Q now mean s fixed colle ction time 
t. In the Poisson limit p — > 0, ( 3.22 ) and ( 3.23 ) reduce to 
the corresponding expressions derived for fixed t in ref. 
p6[ . The expressions for Qe and Qe are obtained by re- 
placing the quantities with suffi x g by the corresponding 
ones with suffix e in (3.21) and ( 3.23D respectively. Note 
that all of these Q-functions are directly proportional to 
the corresponding detector efficiencies for any choice of 
p, thus solving the problem raised, for example, in |2l| ] 
and, for p ^ and fixed t in agreement with [ p6[ . 

We outline in the Appe ndix a meth od of actually eval- 
uating these Q-functions (3.21) and (3.23) for both of e 
and g. 

We have already shown in the last section (at ( p.40[ )) 
that the state of the field observed after the passage of a 
large number A^ of active atoms is the same as the one 
observed after waiting for a long time. We now find the 
atomic Q-functions for a group of a fixed number N —f oo 
of atoms using (3.21) and for a group of atoms passing 
through the cavity in a fixed time i — > oo by evaluating 
( |3.23 ) in the limit t ^- oo. The results of the numerical 
computations are presented in the next section. 

Here we compare the predictions of the two approaches 
in the case of an analytically solvable example. This is 
the example of the microma ser operating without any 
thermal photons (n = in ( ^.2| )) and having an inter- 
action time such that gUnt — ^A/2- In this case the 
Fock state |1) is a trapping state [Q. Hence, if the cavity 
field is initially in the state of vacuum then the state of 
the micromaser is described by two-dimensional matrices 
and this makes the problem analytically tractable. We 
find the Q-functions for the upper and the lower states 
in the limit i — > cxd to be given by the formulae 



-Ve 



-(- 



2p2/32(i _ (^)(c^^ +pl3id^d)d 



and 



where 



l-d + pl3id\{l-d + pl3id)[ai{l -d)+ pM] 
+p[ai{\-d)+p(3id]), (3.24) 



- -T^g / 2p/3i(l-p/3i)(l-rf)d 

'^ \-d + pPid\ 1-d+pPid 

+p/3i(l-d)), (3.25) 

/3i = l-ai =sin2(^), 
A^e. = ^. (3.26) 
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The parameter Nex is the usual so-called pumping pa- 
rameter mentioned in the abstract which gives the aver- 
age number of atoms pumpe d int o the cavity in one cavity 
damping time. Note from ( 3.25 ) that Qg is always neg- 
ative, i.e., sub-Poissonian whereas Qe is definitely sub- 
Poissonian if ai + ppid — rf > 0. For the particular case 
of Poisson pumping statistics the limit p ^ in ( 3.24 ) 
and (3.25) yields 



Qe = 



2VePlN; 



{Nexf3i+ai){Nex(3i + iy 



and 



-2f]e(3iNex 
{NexPl + 1)2 



(3.27) 



(3.28) 



These expressions ( 3.27 ) and ( ^.2^ ) for fixed t ^ oo and 
p ^ are exactly the same as the corresponding ones 
derived in ref. Q . Thus, in the case of Poisson pumping, 
the atoms exiting the cavity in the upper state follow a 
super-Poisson distribution whereas those emerging in the 
lower state have a strictly sub-Poissonian distribution. 

In a similar way the Q-functions in the limit N —f oo 
are found to be given by 



Qe 



-Ve 



2p(ilai{l-dfd 



l-d + pl3id\{l - d + pl3id)[ai{l - d) + pPid] 
+ai(l - d) + pPid\ . (3.29) 



and 



Q9 = 



/2p(3iai{l-d)d 
l-d + pl3id\ 1-d + pPid 



-Vg 



+ Pi{l-d) 



(3.30) 



Note that in this case, i.e., in the case of observations 
made on groups of a fixed large number N of atoms as 
N -^ cxD, both Qe as wel l as Q q are sub- Poissonian. In 
the Poisson limit p ^ 0, ( 3.29| ) and ( 3.3C ) reduce to 



-Ve 



and 



PlNex - 
-Vg 



Wl^lNe: 



1 V [PlNe: 

ai 



l){PlNe: 



OLl) 



t( 



2/3iaiiVe; 



l3iNex^l\{l3iNex + l) 



A 



(3.31) 



(3.32) 



The differences between the atomic statistics in the long 
time and the large number of atoms limits are clearly 
reflected in these various analytical results. 

It is instructive to consider the behavior of the Q- 
functions derived above for the special cases of the micro- 
maser with low and high rates of pumping respectively. 
For the low rate of pumping R (or equivalently high rate k 
of field damping), N^x -^ 0. In this limit (Ng) — > Nrjgfii, 



(iVe) ^ iV77,(l - /3i), Qg -^ -r^gPi, Qe ^ -?7e(l - /3i). 
These values of {Ng) ((iVg)) and Qg {Qe ) correspond 
to a binomial process in w hich the probability of success 



is I3g ([1 — (3g\) (see (4.1) below). Hence it follows that 



the statistics of the atomic populations exiting the cav- 
ity is binomial irrespective of the pumping statistics if 
the observations are made on groups of a fixed number, 
N ^r oo, of atoms. If, on the other hand, the observa- 
tions are made for fixed time, t — > cx) then the desired 
averages are obtained from the corresponding ones for 
groups of fixed number N of atoms by averaging over N 
with the probability that N atoms exit the cavity in time 
t. In this case, for the present example, {Ng) {{Ne)) and 
Qg {Qe) are expected to be Rtrjgfii {Rtr]e[^ — Pi\), and 
—PVaP i (— p ?7e[l — Pi]) [see (4^)]. It is easy to see that 
( |3.24D and (3.25) indeed yield those results in the limit 
Nex -^ 0, SO that d -> 0. 

Next, for Nex —^ oo, that is in the limit of a high 
rate of pumping or equivalently of a low rate of field 
losses (namely for an ideal cavity), one has Qe —^ —1, 
Qe — > — p, Qg —* Qg ^ 0. These are evidently precisely 
the values of the Q-parameters for the atoms entering the 
cavity. This behavior of the Q-parameters can be readily 
understood by noting that in this limit, the steady state 
of the field is the trapping state |1): the atoms therefore 
exit the cavity in their upper states |e) which are the 
states in which they enter the cavity. Hence, in the limit 
Nex — > OO, the statistics of the atoms exiting the cavity 
is exactly the same as the pumping statistics. 

In the next section we investigate in detail the dif- 
ferences between the atomic statistics in the long time 
and in the large N limit by numerical evaluation of not 
only the analytical expressions for the Q-functions for 
the special case oi gtint = t^/V^ at zero temperature dis- 
cussed above, but also for several other gtints defining 
both trapping states and non-trapping states Q at both 
zero and finite temperatures. 



IV. NUMERICAL RESULTS AND DISCUSSION 

In this section we present numerical evaluations of the 
atomic Q functions in the asymptotic limits iV — > oo 
and t —* oo. Since the efficiency factors appear only as 
a multiplicative factor we give the results only for unit 
efficiencies r^e — "Hg = \. 

The parameter Qe {Qg) is obviously the same as Qe 
{Qg) for the case of regular pumping p = 1. We, there- 
fore, examine these parameters first for p = Q (Poisson 
pumping) for this corresponds to a maximum departure 
from the regular input. We display first in Fig. |l| the 
plots of Qe, Qg and Qg, Qg together with Qf as functions 
of the pumping parameter Nex = R/'i.K in the case of a 
Poisson input and for which gtint = i^ / \2 while the sys- 
tem is at zero temperature, n = 0. The exact analytical 
expressions for these Q-parameters were given towards 
the end of the last section. It is clear from the figure 
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that Qe and Qe differ in the entire range of N^x whereas 
Qg and Qg differ significantly for small N^x but approach 
each other as N^x is increased. Notice how it is only Qe 
which closely follows Qf. This is because Qg, Qe and Qg 
all tend to zero with Nex -^ oo while Qe -^ —1. Notice 
that as Nex — > oo either cavity Q —> oo(2k —> 0) or the 
mean entry rate i? — *■ oo (or both) . li R —^ oo the prob- 
ability of there being more than one atom in the cavity 
becomes large and the TC-model Ham iltonian [n6[ must 
replace the JC-Hamiltonian Eq. (2.1). For Nex 



oo 



we must therefore assume the cavity Q becomes large so 
as to maintain any relevance to the actual experiments 
done so far. 

We next examine the effect of finite temperatures on 
the differences in the two kinds of Q-functions by plotting 
in Fig. the same data as in Fig. ^ but for n = 0.1. The 
differences between the two kinds of Q-functions seem to 
decrease with an increase in h. The other point to notice 
is that Qe and Qf were both sub-Poissonian for all Nex 
and were monotonically decreasing to — 1 at zero temper- 
ature as in the Fig. |l|. These behaviors for both Qe and 
Qf and indeed the behaviors of Qg, Qe and Qg evidently 
change dramatically with the increase in temperature 0.2 
K (n = 0.10). We return to this in the Figs. |l^-|l| where 
these effects of temperature are displayed. 

In Fig. H we investigate the effects of the departure 
from the continuous time Poisson statistics on the differ- 
ences between the two limits by plotting the Q-functions 
as functions of Nex for the discrete time binomial dis- 
tribution with p — 0.5 for gtint = 7r/\/2 and n — 0. 
Not too much difference from the plots of Fig. 1 can 
be seen from this significant change of p, except for the 
observation that Qe shows sub-Poissonian behavior for 
the entire range of Nex, whereas for p — 0, it was totally 
super-Poissonian. 

Next we plot in Fig. ^ the Q-functions for another 
value of gtint namely, gtint = 1.54 as a function of Nex 
for n — 0.145 (temperature — 0.5 K) and for the case of 
Poissonian pumping. This is the set of parameters used 
in the experiments of ref. [Effl . Once again we find that 
the two kinds of Q-functions differ considerably and dif- 
fer considerably from Qf for small Nex- It is interesting 
to compare the plot of Qg with the experimental plot, 
Fig. y of ref. |Q. One can check that the theoretical 
plot for Qg given in Fig. here exhibits the minima and 
the maxima at the same positions as the plot of the ex- 
perimental points. However, whereas the second minima 
of the experimental points exhibits sub-Poisson behavior, 
that of the theoretical plot is super-Poissonian and to this 
extent the sub-Poissonian field with Qf = —0.3 observed 
at Nex = 35 in Qg in 1 20 1 may need further scrutiny. 



Note that the numerical values in the two plots, even 
after correcting for the efficiency factor, do not match 
at all. The reason for the discrepancy between the two 
apparently may stem from the fact that there are stray 
fields present in the experiments (cf. ^^) which have not 
been taken into account in the theoretical plot of Fig. 0. 



In Fig. g we plot the same data as in Fig. ^ except that 
h — (zero temperature). In view of the Figs. |l| and H 
where n = and n > the most striking feature of Fig. 
(n — 0) in relation to Fig. |j (n > 0) is that qualitatively 
there is little change, although Qe for Nex ~ 44 is sub- 
Poissonian and approaching Qf from Fig. ^ which is 
also sub-Poissonian for however Nex ~ 35. Note that 
gtint = 1.54 is not a low- lying (in terms of photon number 
n) trapping state. 

Notice how the increase in temperature from zero to 
0.15 K between the Figs. || and ^ pushes all of the Q's 
substantially to the left (lower Nex) so that each of the 
minima of Qe, Qe, Qg, Qg, of which only Qe is sub- 
Poissonian, arise at Nex ~ 44 in Fig. ^ while these are 
pushed to Nex ~ 35 in Fig. H with now Qe just super- 
Poissonian. Apparently the big maximum for the atomic 
Q's for Nex > 42 in Fig. is pushing these minima 
strongly to the left as described. Notice now that this 
effect of temperature for gtint = 1-54 is qualitatively dif- 
ferent from the effect of temperature between Fig. |l| and 
Fig. g. This is apparently due to the fact that in Figs. 
m and 0, gtint = ""/VS which is the one-photon trap- 
ping state condition: the gtint — 1-54 radians does not 
define any trapping state with a small number of pho- 
tons. Thus, with vr == 3.14159 for example, one finds 
that the lowest trapping state has the photon number 
n ^ 4: X 10^°! At zero temperature and Tc = oo (ideal 
cavity) the trapping states are the only equilibrium states 
p3| , plj[ ; but although gtint — 1-54 at zero temperature in 
Fig. g, Tc < oo and equilibrium is expected to lie below 
the lowest trapping state (in terms of photon number) 
as was demonstrated in [|l3|. For lowest trapping states 
with small photon number switching on the temperature 
induces a rather remarkable 'escape' to the larger pho- 
ton numbers which leads to the rather dramatic rises in 
variances evident in the sequence of thresholds for these 
variances shown in the contrasts between the Figs, o and 
shown below for which figures all are taken at low lying 
trapping states. A conclusion remains that there is no 
simple relation between Qf and any of Qg, Qe, Qg and 
Qe at finite temperature while it is only Qe which follows 
Qf relatively closely at zero temperature. 

We next examine, in Figs, o and ^ the differences 
in the two kinds of asymptotic limits as a function of 
gt.nt for Nex = 1 (Fig. |) and Nex = 5 (Fig. § for 
n = 0.1 and for the case of Poisson pumping. These 
figures confirm that these differences depend on Nex for 
any gtint- Thus the curves for Qe, Qg differ from the 
corresponding Qe, Qg curves for all gtint for Nex — 1 
(Fig. @), while their differences decrease as the values of 
Nex is increased to 5 (Fig. ||). The Fig. g shows these 
differences for the experimentally realized value of Nex = 
30. In the Fig. |9| we repeat the data of Fig. at n = 
(zero temperature). Low- lying trapping states arise 
at gtint = tt/v^, 7I'/\/3, 7r/2, n/^/b, 7r/-\/6 for example 
for r = 1 in the trapping state condition gtint ^/n + 1 = 
TTT. For r — 1,2,3 etc. they arise for gtint = 7r/2, 7r/4, 
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7r/6 and so on, for example. Still Fig. shows little 
qualitative change from Fig. || (n = 0). 

For the Fig. 00^ we note that the value of gtint , gtint = 
7r/V2, in the Fig. |l| is a small value for a trapping state. 
In the Fig. nfl we choose n — (zero temperature) and 
gtint = 7r/2 radians. This is a trapping state with photon 
number n = 3 also a small value. The behaviors in Fig. 

10 are qualitatively the same as in the Fig. ^ In the Fig. 

11 the temperature is increased so that h = 0.1. The 
qualitative change compared with Fig. |l^ follows the 
qualitative change between Fig. n (zero temperature) 
and Fig. || (finite temperature). 

For the Fig. [l^ we choose gtint — 7r/-\/l9 so that the 
trapping state photon number is increased to n = 18, 
larger than in Figs. ^ and and fol and |ri|. The Fig. 
O shows the same data for n — 0.1; the change here is 
the emergence of the threshold for Qe, Qg just beyond 
A^ex = 60. 

The Figs. |lj (a) and (b) plot the dependence of Qg 
and Qe respectively on the collection time t for gtint = 
7t/\/2. The chosen t values are 1, 5, 10, 100, 1000 and oo 
times the cavity damping time and the plots are given as 
functions of N^x- These results compare with the Fig. H 
in |4^ evaluated for gtint — 1-54 and plotted as a function 
of collection time t. Finally the Figs. 11^ (a) and (b) 
plot the dependence of Qg and Qe respectively on N the 
fixed number of atoms collected as a function of N for 
gtint = 1-54 and n = 0.1. The chosen N values are 20, 
50, 100, 500, 1000, and oo. 

The reason for the differences in the atomic Q- 
functions in the limit of a large number of atoms and that 
in the limit of long time for N^x <C 1 can be easily seen. 
These limits correspond to high cavity damping or a low 
rate of injection of the atoms. In either case the cavity 
field decays to a state of thermal equilibrium during the 
time interval between the arrival of two successive atoms 
into the cavity. Hence each atom sees the same field, i.e. 
the thermal field, on its entry into the cavity. The prob- 
ability of an atom exiting the cavity in a given state is, 
therefore, the same for every atom. In other words, there 
is no correlation between the states of two atoms exiting 
from the cavity and so the probability that an atom exits 
in a given state follows a binomial distribution. Let (3g be 
the probability that an atom exits from the cavity in its 
ground state. If one observes groups of a fixed number 
N of atoms exiting from the cavity then the number iV^ 
of atoms detected in the ground state is governed by a 
binomial distribution for which the probability of success 
is (3g. Hence 

{Ng) = N(3g, 



{NgiNg~l))=NiN- 



iK 



(4.1) 



As a consequence of (4.1) it follows that Qg = —fig- 
The same relations hold for the Q-function Qe for the 
upper state with f]g replaced by 1 — /3g. For the ex- 
actly solvable case of the micromaser discussed above, 
/9g = sin^(7r/-\/2) « 0.36. Hence Qg and Qe in this case 
are predicted to be approximately —0.64 and —0.36 re- 



spectively for Nex close to zero. This is in agreement 
with the plots presented in Figs. |l|. 

If, on the other hand, one observes the groups of atoms 
for a fixed time interval t then the number N of atoms 
entering the cavity in that interval is governed by the 
pumping statistics. The expressions corresp ondi ng to 
(4.1) in this case are obtained by averaging (4.1) over 
the number TV of atoms. Now, the pumping statistics 
has been assumed to be a binomial distribution for which 
the probability of success is p. Hence (N) = Rt and 
{N'^) - (7V)2 = Rt{Rt - p). On averaging {^) over N 
we get 



(Ng) = RtP, 



3' 



{Ng{Ng^l))=Rt{Rt-p)l3l (4.2) 



As a result, Qg = —pf3g. The same arguments hold for 
Qe with (3g replaced hy 1 — Pg. In particular, for the 
Poisson pumping p — 0, it is expected that Qg = Qe — 
for Nex close to zero. This is clearly in agreement with 
the plots presented above. 

As Nex increases, the correlations between the succes- 
sive atoms start dominating. In particular for Nex -^ oo, 
i.e., in the limit of an ideal cavity, the steady state of the 
cavity field approaches a trapping state and all the atoms 
exit in the excited state which is the state in which they 
enter the cavity. Hence, in this case, it is expected that 
Qe ^ —I, Qe ^ —p, Qg ^ Qg ^ 0. With the exception 
of Qe whose asymptotic, large Nex behavior could not 
actually be shown in the figures this behavior is borne 
out by the figures presented above. 
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APPENDIX A: 

In this Appendix we outline a met hod f or evaluating 
the Q-functions given by ( 3.21 ) and ( 3.22 ). Note that 
those two Q-parameters can be written as 



2% 



Trf[FgPs 
1 



rTr 



[l-AB) 



'V K{\-AB)) 
AFgPss -rigKTrf[FgPs 



(Al) 



where K ^ N , A ^ Ap, B = Fa ioi Q = Qg and K = Rt, 
A = Ap,B = Fo for Q = Qg. 

To evaluate (Al), it is convenient to work in the eigen- 
representation of the operator AB. Let {pi} and {pi} be 
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the right and left eigenfunctions of AB corresponding to 
the eigenvalues {Ai}, i.e., 



ABpi ^ XtPt, PiAB = Xipi 



(A2) 



We assume that the eigenvalue Aq = 1, i.e., po is the 
steady state density matrix of the micromaser (p^s = Po)- 
On taking the trace on both sides of the first of the 
equations in ( A2) and on using the trace conserving prop- 
erty of A and i?, i.e., the relations Tr{Ap) = Tr{p), 
Tr{Bp) = Tr{p), which are a consequence of the iden- 
tification of A and B with the operators Ap, Ap, Fq, Fq 
and the trace conserving property ( 3.19 ), it follows that 
Tr(pi) = XiTr(pi), and since Aq is the only eigenvalue 
assumed to be unity, this leads to the result 



Tr{p,)^0, if 1^0. 



(A3) 



Now, the eigenstates {pi} and {pi} are orthogonal with 
respect to the trace as a scalar product. Let us assume 
that they are orthonormalized so that 



Tr 



yPtPj) 



Si 



As a consequence of (|A4D, we can write 



where 



AB{pi)^J2^fM 



C!^^Tr(p,ABp, 



(A4) 



(A5) 



Now we can express the averages and the Q-functions in 
terms of the eigenstates of AB using ( A6). 

Since Tr{AFgp) = Tr{Fgp), it follows on using (^) 
that 

Tr{Fgp) = Co^Q. (A7) 



In the same way, it can be shown that the expression 
(Al) for Q reduces to 



-W9<-oo 



2^3 

C^nnK ^ L( 



r<9 r<9 



i:[7r&("-W'--' 



-\l-p(l-\i)] 



K 



(A8) 



The Q-functions can thus be computed by evaluating the 
right and the left eigenstates and eigenvalues of the op- 
erator AB, where A and B are defined after (Al). 
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FIG. 1. The Q-parameters, Qg, Qe, Qg, Qe, and Qf as 
functions of N^x for gtint = 7r/\/2, n = and for Poissonian 
pumping of the excited atoms into the cavity {p = 0). Qf 
is multiplied by 0.1 in order to compare it with the atomic 
Q-parameters. 




FIG. 2. The same as Fig. hi, but for non-zero thermal pho- 
tons (n — 0.1). 
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FIG. 3. The same as Fig. hi, but for binomial pumping of 
the excited atoms into the cavity with excitation probabihty, 
p = 0.5 (deviation from Poissonian pumping). 




FIG. 5. The same as Fig. |, but for n = 0. 




FIG. 4. The Q-parameters, Qg, Qe, Qg, Qe, and Qf as 
functions of Nex for gti„t = 1.54 and n = 0.145 and for Pois- 
son pumping of the excited atoms into the cavity {p = 0). 
These parameters correspond to the Fig. of Mt . 
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FIG. 6. The Q-parameters, Qg, Qe, Qg, Qe, and Qf as 
functions of gtint when n — 0.1 and N^^ = 1 for the case of 
Poissonian pumping of the excited atoms. 
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FIG. 7. The Q-parameters, Qg, Qe, Qg, Qe, and Qf as 
functions of gtint when n — 0.1 and Nex ~ 5 for the case of 
Poissonian pumping of the excited atoms. 
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FIG. 9. The same as Fig. ^ but with n = 0. 




FIG. 8. The same as Figs. | and |, but for Nex 
Computational difficulties have prevented us following 
17 < gU^t < 30. 
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10. The same as Fig. hi but for gtint = ir/y/A = 7r/2 
= 0. 
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FIG. 11. The same as Fig. hd but for finite temperature 
(n = 0.1). 
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FIG. 13. The same as Fig. h2 but for finite temperature 
(n = 0.1). 
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FIG. 12. The same as Fig. ^ but for gti„t = 7r/\/T9 and 
n = 0. 
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FIG. 14. The dependence of (a) Qg and (b) Qe on the 
collection time t for gtint ~ 1-54 and n = 0.1. The t values 

are 1 ( ), 5 (• ■ ■), 10 ( ), 100 ( ), 1000 ( ), 

and cx) ( — ). 
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FIG. 15. The dependence of (a) Qg and (b) Q^ on the 
number of atoms collected, A*' for gtint = 1.54 and n — 0.1. 

The iV values are 20 ( ), 50 (•••), 100 ( ), 500 ( ), 

1000 ( ), and c» (— ). 
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